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ABSTRACT 

The stellar initial mass function (IMF) plays a crucial role in determining the 
number of surviving stars in galaxies, the chemical composition of the interstellar 
medium, and the distribution of light in galaxies. A key unsolved question is whether 
the IMF is universal in time and space. Here we use state-of-the-art results of stel- 
lar evolution to show that the IMF of our Galaxy made a transition from an IMF 
dominated by massive stars to the present-day IMF at an early phase of the Galaxy 
formation. Updated results from stellar evolution in a wide range of metallicities have 
been implemented in a binary population synthesis code, and compared with the ob- 
servations of carbon-enhanced metal-poor (CEMP) stars in our Galaxy. We find that 
applying the present-day IMF to Galactic halo stars causes serious contradictions with 
four observable quantities connected with the evolution of AGB stars. Furthermore, 
a comparison between our calculations and the observations of CEMP stars may help 
us to constrain the transition metallicity for the IMF which we tentatively set at 
[Fe/H] « — 2. A novelty of the current study is the inclusion of mass loss suppres- 
sion in intermediate-mass AGB stars at low-metallicity. This significantly reduces the 
overproduction of nitrogen-enhanced stars that was a major problem in using the 
high-mass star dominated IMF in previous studies. Our results also demonstrate that 
the use of the present day IMF for all time in chemical evolution models results in the 
overproduction of Type 1.5 supernovae. More data on stellar abundances will help to 
understand how the IMF has changed and what caused such a transition. 
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1 INTRODUCTION 

The IMF is one of the most important factors influenc- 
ing the evolution of galaxies. In particular, possible tran- 
sitions of the stellar IMF from that of the first (metal-free) 
stars to that of the present-day (metal-rich) stars have been 
given much attention. The typical first generation stars are 
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thought t o be predominantly of high-mass (e.g. ~ 40Mq ac- 
cording to lHosokawa et al.ll2011h from studies of star forma- 
tion in metal-free environments. However, the actual form of 
the IMF for such stars is poorly un derstood and is still con- 
troversial (see e.g. iGreif et al.ll2012l . and references therein). 
In addition, there are few studies dealing with the transi- 
tion of the IMF during the early epoch of the Galaxy. If the 
IMF is not universal then a question arises - how and when 
did the IMF change to the present- day IMF which is biased 
towards low-mass stars (< Q.8Mq) (Sal peteJl 19551 )? 
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Extremely metal-poor stars (EMP stars) observed in 
the Galaxy can help to answer this question because we 
can estimate the number of former AGB stars that con- 
tributed to the change of the surface ab undances evident in 
these stars. Spectroscopy of EMP stars dBeers fc Christliebl 
2005) reveals that carbon-enhanced stars are very com- 
mon amongst EMP sta rs (20 % or more l|Suda et al.l 
l201ll ; ICarollo et al.l [20121 )) when compared to the higher- 
meta llicity subgiant CH stars (~ 1 %) (|Luck fc Bond! 
Il99ll ). The large frequency of carbon-enhanced metal- 
poor (CEMP) stars is argued to be a consequence of 
binary mass transfer from AGB stars th at have under- 
gone carbon dredge-up into their envelopes llFuumoto et al 



2000; Suda et al. 2004; iLucatello et al.ll2005bl ; lKomiva et al 



20071 ). The mechanism for carbon enhancement, the third 
dredge- up (TDU), is the currently accepted theory that ex- 
plai ns the origin of classical Ba II stars and CH stars (see, 
e.g. lMcClure fc WoodswOTthTl 990). In the same manner, the 
origins of some CEMP stars are thought to be explained 
by the binary scenario, i.e., the observed stars have binary 
companions that are now unseen white dwarfs. One of the 
advantages of this scenario is that we do not expect all the 
observed EMP stars to show carbon enhancement on their 
surface. This is qualitatively consistent with the recent dis- 
covery of EMP stars in the lowest metallicity range with- 
out large enhancements of carbon (Caf fau et aj]|201ll ). Bi- 
narity has been statistically confirmed for the population 
of CE MP stars that show enhancements of s-p rocess ele- 
ments (|Suda et alj |2004| ; iLucatello et al.1 l2005bl ). as is ex- 
pected from theory. These objects are called CEMP-s stars, 
in contrast to the other subclass of CEMP stars, the CEMP- 
no stars, which do not show any enhancement of s-process 
elements. The origin of CEMP- no stars is still controver- 
sial (see, e.g., iNorris et al.ll2013l ). Observationally, there is 
another class of EMP stars that show enhancements of ni- 
trogen greater than the enhancement of carbon. These stars 
are called nitrogen-enhanced metal-poor (NEMP) stars and 
are thought to be produced via binary mass transfer from 
intermediate-mass AGB stars which have converted C to N 
via hot bottom burning (HBB) during the AGB phase. 

Previous works exploring the p ossibility that the IMF 
was different in th e early universe (|Lucatello et al.l [2005a; 
IKomiva et al.ll2007l) dealt o nly with the statistics of CEMP 
stars. ( Komiva et al.l l|2009l ) also pointed out that the differ- 
ent IMF is suggested by the number of low-mass star sur- 
vivors estimated from the survey areas of HK and Hamburg- 
ESO surveys.) In particular, t he IMF which was pr oposed to 
have a peak at around WM@ |Komiva et al.ll2007l ) has been 
criticized due to th e possible overpro duction of NEMP stars 
|lzzard et al.ll2009l ; IPols et al.ll2012l ). However recent work 
on pulsation driven mas s loss in AG B stars at metallici- 
ties below [Fe/H] ~ -2.5 (|Woodll201ll ) suggests that NEMP 
star formation should have been suppressed in the early Uni- 
verse. That study showed that the growth of pulsation dur- 
ing the AGB phase is too small to trigger dust-driven mass 
loss in low-metallicity, intermediate-mass stars. With mass- 
loss suppressed, the degenerate cores of these stars would 
have grown to be massive enough for the stars to explode as 
type 1.5 supernovae, i.e., the mass of their cores wou ld have 
grown to approach the Cha ndrasekhar mass limit l|Arnettl 
Il969l ; llben fc Renzuiilll983D . The purpose of this study is 
to explore the possibility that the IMF is different at low- 



metallicity by including as many ingredients related to stel- 
lar evolution and binary evolution as possible. In particular, 
this is the first study to in clude characteristic A GB yields 
at [Fe/H] < -2.5 (see, e.g.. iFuiimoto et aljlioool ). 



2 BINARY POPULATION SYNTHESIS OF 
EMP STARS: CODE DETAILS 

We use the Monte Carlo method to simulate the evolu- 
tion of binary systems by generating binaries with a range 
of initial primary mass, secondary mass, and binary pe- 
riod. The mass of the primary is described by a log- 
normal function or power-law function. The former is de- 
fined by dN/dlogm oc exp(-(log(m/^)) 2 /(2cr 2 )). The me- 
dian mass and the dispersion are set at (/i (in M@), a) — 
(5, 0.6), (10, 0.4), (20, 0.45), (30, 0.5), and (50, 0.6), re- 
spectivel y, in order to reprod uce CEMP-no/CEMP-s ratio 
following IKomiva et al.l {2007). Power-law functions are de- 
fined by dN/dlogm tx m~ x , where we choose x = 1.35, 
0.85, 0.35, and 0.0. The lower and upper cut-off mass 
are set at O.O8M0 and 2OOM0, respectively. The mass 
of the secondary is subje ct to the mass ratio funct ion 
which was taken from the iDuquennov fc Mayor! (ll99lT) or 



iRaghavan et all (|2010l ). For the lRaghavan et aLp|2010l ) dis- 



tribution we adopted a flat distribution by ignoring the 
preference for like-mass pairs. For the period distribution 
of binaries, we adopted a fitting formula f or nearby main- 
sequence stars |Duquennov fc Mavorlll99ll ; IRaghavan et al.l 
|2010| ; lRastegaevll2010l )~ 

The simulations generate binary systems for metallic- 
ities in the range —6 ^ [Fe/H] —1 in steps of 0.5 dex 
in each simulation. The total number of binary systems for 
a given metallicity is typically one million, which produces 
more than 1,000 observable giants. We note here that the 
simulations are compared to observations of giant stars (Fig- 
ure [TJ. Giants were chosen so that the uncertainty concern- 
ing the depth of the convective envelopes for main-sequence 
stars is reduced, and because the uncertainty in the ele- 
mental abundances derived for dwarfs is considerable. We 
assume that the stars with M ^ 0.83Mq evolve to the final 
stage of their evolution. 

In our simulations, we consider, for the first time, 
another channel of CEMP star formation, i.e., hydrogen 
ingestion into the helium flash convective zone for AGB 
stars at [Fe /H] ^ -2.5 (He-Fla s h Dri ven Deep Mixing; 
He-FDDM) (|Fuiimoto et al.1 Il990l . |2000| ) . Previous studies 
have not taken t his channel into account (jlzzard et al. I l2009l ; 
IPols et al.|[2012l ). Surface chemical abundances of primary 
stars are taken from stellar models after the He-FDDM 
at the beginning of the therm ally pulsing AGB (TP- AGB) 
phase |Suda fc Fuiimotoll2010l) . or after the TDU during the 
TP- AGB phase. 

The surface abundance evolution for all the stars was 



taken from the literature ( 


Suda & Fuiimotoll2010l; iKarakas! 


l20ld; iGil-Pons & Dohertvl 


20ld; IGil-Pons et al.1 20131). For 



models undergoing the He-FDDM event, we took the chem- 
ical abundances just after the dredge-up since this event 
dominates the change of the surface composition. For the 
TDU event, we took t he avera g e abu ndances of chemi- 
cal yields provided by IKarakas] (|2010l ). For models with 
[Fe/H] < —2.5 with the mass range responsible for the TDU, 
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we adopted the same yields as in iKarakasI (|201Ch by adopt- 
ing the values at [Fe/H] = —2.3 of the same mass. We also 
considered the effect of carbon enhancement by super- AGB 
stars that host oxygen-neon-mag n esium cores dRitossa et al.l 
1 19991 ; iGil-Pons fc Dohertvl l201ol ; iGil-Pons et al l l2013h . al- 
though it turned out to be a minor contribution to CEMP 
stars in our simulations. A linear interpolation is made for 
the surface abundances between the model grids of initial 
mass and metallicity. 

Binary mass transfer events are treated as instanta- 
neous episodes. They are described by the masses of the 
binary components and the binary separation assuming Ke- 
pler rotation without any orbital ellipticity. If the binary 
separation is larger than the sum of the Roche radius of 
the primary star and the radius of the secondary star dur- 
ing the main sequence phase, wind accretion is assumed to 
be at play. Otherwise, Roche lobe overflow (RLOF) applies 
to the system. In this study, the occurrence of mass loss in 
AGB stars is controlled by a parameter introduced to take 
into accoun t of the poss ible suppression of mass loss at low 
metallicity (|Woodl l201ll l . Here we assume that stars with 
5 < M/Mq < 8 and [Fe/H] < —2.5 do not eject any envelope 
and explode as Type 1.5 supernovae, which is chosen to be 
consistent with the measur ed abundances using the SAGA 
database jSuda et al. 2008). The SAGA database includes 
451 giant stars with carbon abundances derived from spec- 
tra with a resolution of R — A/AA > 20000. The database 
covers the data published up to 2011 (August, 2012 update). 

For wind accretion episodes we assume that the mass 
accretion onto the secondary star is subject to Bondi-Hoyle 
accretion with the constant stellar wind velocity of 20 km/s 
from the primary star. The amount of mass lost from the 
primary is estimat ed from initial-m ass final-mass relations 
for white dwarfs (|Han et al. These treatments are 



the same as in our previous work (IKomiva et al.ll2007) and 



the details are given in the reference l|Komiva et alj 



llzzard et al]|2009l ; lEggletorj[l983l ; iHurlev et alj|2002h 



20071 ; 



For RLOF, the Roche radius is calculated from a com- 
monly used approximation , Ri/a = 0.49q 2//3 /(0.6g 2 / 3 + 
ln(l + g 1//3 )), where Ri, a, and q = M1/M2 are the radius of 
a primary, bina ry separation, and mass ratio, respectively, 
teggletorJ ll983h.' If the condition for RLOF accretion is sat- 
isfied, we assume a constant m ass accretion of 0.0 5Mp) onto 
the secondary occurs, following llzzard et al.l (2009). There is 
however no supporting evidence for this exact amount. An 
important point here is that the common envelope phase 
can be avoided only for binaries that have a mass ratio very 
close to unity. Thus the effect of RLOF accretion is actually 
always negligible in this study, and most binaries undergo- 
ing RLOF are excluded from the observable population due 
to merging. 

The definition of CEMP giants is given by [C/Fe] > 0.7 
l|Suda et al.l l201ll . see also Aoki et al. 2007) in the mass 
range 0.82 to O.83M0, considering the stellar lifetimes in the 
Galactic halo population (comparable to the age of the Uni- 
verse). Self-pollution by AGB stars is excluded due to their 
negligible contribution to the population compared with red 
giants. We designate CEMP star candidates as NEMP stars 
if the initial mass is in the range in which HBB operates. 
This mass range is a parameter of our model and is set at 
4.5 < M/Mq < 8 in our fiducial models. For the defini- 
tion of CEMP-5 and CEMP-no stars, we use parameters for 



mass and metallicity. We tentatively define a critical mass, 
set at 3.5M0, below (above) which the model stars become 
CEMP-s (no) stars during the AGB phase. For CEMP- 
no stars, it is required that the model stars should have 
[Fe/H] ^ —2.5 to be consistent with observations. Note that 
the origin of CEM P-no stars is still controversial (see, e.g. 
iNorris et al.ll2013T ) and will be discussed i n a separate paper. 
See also discussions in ISuda et all l|2004l ) for more detail. 



3 RESULTS AND DISCUSSION 

The frequency of CEMP stars as a function of metallicity is 
given in Figure [1] for observations and for our models with 
different adopted IMFs. The CEMP fraction CEMP/EMP 
and the fractions CEMP-no/CEMP and NEMP/CEMP 
from our models are shown in Table[T]for a selected set of pa- 
rameters. These ratios are averaged values from the SAGA 
database for metallicities below [Fe/H] = —2. They are ob- 
tained by taking a weighted average of the ratios in each of 
the metallicity bins below [Fe/H] = —2, with the weight be- 
ing the number of stars in each bin. Model B is our fiducial 
model as it satisfactorily reproduces the observations. It was 
chosen to reproduce the distribution of carbon abundance as 
shown in Figure(2] It is to be noted that the CEMP fractions 
of stars taken from the SAGA database can be affected by 
selection bias, although they are consiste nt with the values 
derived for other h omogeneous datasets |Yong et al.ll2013l ; 
ICarollo et ai1l2012T l. 

The results presented in the top panel of Fig. Q] demon- 
strate that the CEMP fraction for EMP stars is well repro- 
duced by the high-mass star dominated IMF (model B), 
while the Salpeter IMF (model G) produces a very low 
CEMP fraction compared w ith observations. This is con- 
sisten t with previous work |lzzard et al.l |2009| ; I Pols et ah! 
2012), although we note that those studies did not consider 
the effects of stellar evolution at [Fe/H] < —2.5 and as- 
sumed that all the stars were born as binaries. The CEMP 
fraction produced by the model with the flat IMF (model J) 
is also below the observed value for [Fe/H] ^ —2. These re- 
sults show that in order to reproduce the CEMP fraction at 
[Fe/H] ^ —2, the EMP population should be dominated by 
massive stars that end their lives as type II supernovae i.e. a 
high-mass star dominated IMF is required. For more metal- 
rich populations ([Fe/H] > —2) such as the inner halo pop- 
ulation, a flat IMF or low-mass star dominated IMF gives 
a better agreement with observations. This implies that an 
IMF transition occurred at [Fe/H] ~ —2, perhaps in associa- 
tion wiHi_y2£^trwtirreformation process of our Galaxy (see 
also lYamada et al"1l2013I ). In the bottom panel of Fig. [1] we 
present a model with a transition of the IMF at [Fe/H] = — 2. 

Here we briefly mention the parameter dependencies of 
the simulations, although these are thoroughly discussed in 
a separate paper (Suda et al., in prep.). For the choice of 
the period di s tribut ion function, the distribution derived 
by iRastegaevI (|201(t) gives preference to short period bi- 
naries, which results in efficient channels for the common 
envelope phase and CEMP star formation. We adopted 
this distribution because it is the only study that consid- 
ers a population of metal-poor b inary stars. Using the other 
two binary period dist ributions l|Duquennov fc Mavorlll99ll ; 
Ragh avail et al.ll2010t ) will decrease the fraction of CEMP 
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stars by up to ~10% for the high-mass IMF and for 
the low-mass IMF. Our fiducial model adopts other input 
physics and parameters so that the number of CEMP stars 
can be maximized with reasonable choices. The choice of 
the mass ratio function does not s ignificantly alt e r the frac- 
tion of CEMP sta rs as discussed in lKomiva et al.l (|2009t ) and 
iPols et all (|2012h . The boundary masses for the occurrence 
of HBB and suppression of mass loss significantly affect the 
fraction of NEMP stars. These values are still open ques- 
tions and cannot b e determined at this stage fsee lPols et al.l 
l2012l ; IWoodll201ll ). The mass boundary for the progenitors 
of CEMP-no stars is also an open question and can change 
the ratio CEMP-no/CEMP. 

There are four observable constraints on the simula- 
tions. Three of them are shown in Table Q3 and the other 
one is the [a/Fe] abundances of EMP stars. Firstly, it can 
be seen that the CEMP fraction in our models can never 
reach as high as the observed 25 per cent when using a 
low-mass IMF (see also Fig. [TJ. In order to check the ro- 
bustness of the conclusion, we have tested extreme assump- 
tions for the model parameters, which include (1) a binary 
rate of 100 %, (2) excluding the contribution of supernova 
binaries because the speed of companion stars could have 
been higher than the velocity dispersion of the host ha lo af- 
ter the disruption of binaries l|Tauris fc Takensl ll998). and 
(3) the carbon abundance of primary star enhanced by a 
factor of 10. The superposition of all these assumptions 
still fails to explain the observed CEMP frequencies at low- 
metallicity when a low-mass IMF is used. A significant in- 
crease in the CEMP fraction can only be obtained by as- 
suming that all the EMP stars are born as binary systems, 
although this gives the CEMP fraction well below 20 % at 
[Fe/H] = —2.5. A 100% binary fraction is however highly 
unlikely dpuquennov fc Mavorlfl99ll ; iRaghavan et al.ll2010l ; 
lRastegaevll2010h . 

Secondly, we find that the ratio CEMP-no/CEMP is 
extremely low for the models using the low-mass IMF. To 
obtain a ratio CEMP-no/CEMP = 0.5, we have to set the 
lower mass boundary for the formation of CEMP-no stars 
as small as ~ IMq. However, this then leads to a mismatch 
with the observed ratio CEMP-s /EMP. Thus it is evident 
that other sources of CEMP-no stars, for example stars that 
formed out of an ISM pre-enr iched with carbon, are needed 
if we adopt a low-mass IMF l|Mevnet et al.ll2010l ). We note 
however that any source of pre-enrichment would require 
more massive stars, i.e., a massive-star dominated IMF, al- 
though its shape depends on the efficiency of low-mass star 
formation out of ejecta from the massive stars. 

Thirdly, the ratios of NEMP to CEMP stars in models 
with low-mass IMFs are too small to be compatible with the 
observations. The o bserved ratio is around ~ 0.1 in t he cur- 
rent available data ijSuda et al.ll201ll : |Pols et al.ll2012ri. while 
it wa s previously thought to be much smaller |johnson et alj 
2007). The mass range of NEMP star progenitors would 
need to be 2.5 ^ M/M@ ^ 8, without the assumption 
of mass loss suppression, to obtain a NEMP/CEMP ratio 
~0.1. This mass range would need a higher efficiency of 
HB B at low metallicity as has be en reported by some stud- 
ies (|Campbell fc LattanziolfeoOSf ). However, this range will 
also reduce the CEMP fraction dramatically, and hence can- 
not coexist with the requirement for the large frequency of 
CEMP stars. 



Finally, the number of the progenitors of type 1.5 super- 
novae dominates over, or is comparable to, that of type II 
supernovae when a low-mass IMF is used, which is not sup- 
ported observationally since Galactic halo stars do not show 
the low values of [q/Fe] that would be expected from type 
1.5 supernovae (|Nomoto et al.lll984l ). Thus our simulations 
arrive at the robust conclusion that population simulations 
using low-mass IMFs cannot match the observational con- 
straints given by Galactic EMP stars. 

Our results add to previous evidence for a non- 
universality of the IMF. Other evidence comes from studies 
such as the population of white dwarfs (|Adams fc Laughlinl 
19961). those involving simulations of star formation (see, e.g. 
Larson! l2005h. an d obs ervations of early- type galaxies (see, 
e.g. van Dokkuml l2008'). Our study shows that the observed 
population in the Galactic halo can be explained by a transi- 
tion of the IMF from one biased towards high masses to one 
biased toward lower masses at a metallicity of [Fe/H] ~ —2. 
We propose that the statistics of CEMP and NEMP stars 
in the entire range of metallicity in our Galaxy and dwarf 
galaxies in the local group should be investigated to reveal 
the star formation history of galaxies. 



4 CONCLUSIONS 

We have compared the results of our binary population syn- 
thesis model with the observations of metal-poor halo stars 
in the Galaxy. We find that a low-mass biased IMF, similar 
to the present one, is incompatible with the four observable 
quantities related to the CEMP subclasses. On the other 
hand, considering a high-mass biased IMF, together with the 
occurrence of SN 1.5 in a certain mass and metallicity range 
(supported by recent simulations of pulsation-driven mass- 
loss in low metallicity stars), a more reasonable agreement 
between models and observations is found. Our results imply 
that the transition of the IMF occurred around a metallicity 
of [Fe/H] ~ —2, which should be tested by star formation 
theory and observations. These results also give important 
clues in the most metal poor regimes as to the mass and 
metallicity limits of stars that experience HBB, TDU and 
stellar winds, or that allow the formation of CEMP (-s and 
-no) stars, NEMP stars, SN 1.5 and SNII. Our understand- 
ing of these limits has been traditionally hampered by the 
scarcity of observations. We hope these results can help to 
improve the understanding of the chemical evolution of our 
Galaxy. 
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Table 1. Characteristics of simulated EMP stars with 
the different IMF 



ID 


(/i, tr) or x 


CEMP d 
EMP 


CEMP-no d 
CEMP 


NEMP d 
OEMP 


A 


5, 0.6 


0.19 


0.10 


0.036 


B" 


10, 0.4 


0.28 


0.24 


0.089 


C 


20, 0.45 


0.22 


0.29 


0.092 


D 


30, 0.5 


0.19 


0.29 


0.104 


E 


50, 0.6 


0.19 


0.28 


0.082 


F b 


0.79, 0.51 


0.037 


0.00 


0.072 


G c 


1.35 


0.08 


0.02 


0.010 


H 


0.85 


0.10 


0.02 


0.013 


1 


0.35 


0.12 


0.04 


0.016 


.1 


0.0 


0.14 


0.07 


0.032 


Observed 




0.25 


0.5 


~0.1 



a iKomiva et al.l (120071 1 



*> [ Lucatello et al.l j2005af) 
c ISalpeterl lll955[) 

d Averaged values considering the metallicity distribu- 
tion function of the metal-poor stars for [Fc/H] —2. 
See text. 
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Figure 1. Top panel: the CEMP fraction CEMP/EMP for models with a log-normal IMF with the peak mass at WMq (filled triangles, 
red), a fiat IMF (filled squares, black), and a Salpeter IMF (filled circles, b lue). The model n ames in the top-right corner correspond to 
those in Table [T] The observed CEM P fraction from th e SAGA database dSuda et al.| [2008) is shown by filled inverted triangles (red) 
and the observed CEMP fraction from lYong et al" I d2013h is sh own by open square s (orange). The CEMP fractions for the outer halo and 
the inner halo, kinematically selected from the SDSS sample (Carollo et al. 2013), are shown by filled and open diamonds, respectively 
(both green). The error bars for models and the SAGA sample are based on the bootstrap method with 1000 trials. For [Fe/H] ^ —4 
(shaded area), the number o f stars in the obse rved sample is less than 10 in each bin and statistically unreliable. The criterion for CEMP 
stars is set at [C/Fe] ^ 0.7 llSuda et alj|201lh for both models and observations. Bottom panel: the CEMP fractions for a model where 
the transition of the IMF at [Fc/H] = —2 is taken into account by combining the results for the high-mass IMF (model B) and the 
low-mass IMF (model G). The combined model is compared with the same set of observations as in the top panel. The hatched area 
around [Fe/H] = —2 represents the possible transition metallicity of the IMF from the high-mass star dominated IMF to the present-day 
IMF. 
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Figure 2. The comparison of the observed carbon abundance distribution and the distribution of carbon abundance for the fiducial 
model B. Histograms show the observed distributions for CEMP giants taken from the SAGA database, as used in Fig. [T] Model results 
are shown by solid lines where the number of stars is arbitrarily scaled. 



